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Abstract 

We explore the quantum-classical crossover of two coupled, identical, superconducting quantum inter- 
ference device (SQUID) rings. The motivation for this work is based on a series of recent papers. In 
[1] we showed that the entanglement characteristics of chaotic and periodic (entrained) solutions of the 
Duffing oscillator differed significantly and that in the classical limit entanglement was preserved only 
in the chaotic-like solutions. However, Duffing oscillators are a highly idealised toy system. Motivated 
by a wish to explore more experimentally realisable systems we extended our work in OlS] to an 
analysis of SQUID rings. In [5] we showed that the two systems share a common feature. That 
is, when the SQUID ring's trajectories appear to follow (semi) classical orbits entanglement persists. 
Our analysis in [3] was restricted to the quantum state diffusion unravelling of the master equation - 
representing unit efficiency heterodyne detection (or ambi-quadraturc homodync detection). Here we 
show that very similar behaviour occurs using the quantum jumps unravelling of the master equation. 
Quantum jumps represents a discontinuous photon counting measurement process. Hence, the results 
presented here imply that such persistent entanglement is independent of measurement process and 
that our results may well be quite general in nature. 

1 Introduction 

In this work we extend the results of a recent paper [3] where we investigated the entanglement properties 
associated with the quantum classical crossover of two coupled superconducting quantum interference 
device (SQUID) rings (comprising of a thick ring enclosing a Josephson junction). Here we present a 
small but significant extension of the series of papers [IH2] which forms a small part of a much larger 
body of of work - or example see [4l-[TT]). In order to avoid too much repetition of text please see [iHSllS] 
and references therein for a more detailed introduction to the subject. Here we present a brief summary 
of [3] and our result. 

In [3] we demonstrated that two coupled SQUID ring's can exhibit entanglement that persists even 
in the correspondence limit. In order to obtain these trajectories we used the quantum state diffusion 
unravelling of the master equation and followed a strategy that has seen a lot of success with classically 
chaotic systems [5]. However - there are an infinite number of ways to unravel the master equation. Hence, 
a natural concern that arises is that this result might be unravelling dependent. Here we show that very 
similar behaviour occurs using the quantum jumps unravelling of the master equation. Quantum jumps 
represents a discontinuous photon counting measurement process. 
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Here our interest lay in understanding how the quantum mechanical phenomena of entanglement 
would change as the coupled system approached the classical limit. We showed "that the entanglement 
characteristics of two 'classical' states (chaotic and periodic solutions) differ significantly in the classical 
limit. In particular, we show[edj that significant levels of entanglement are preserved only in the chaotic- 
like solutions" [1]. In [3] we extended this investigation to study the entanglement characteristics in the 
quantum-classical crossover of two identical coupled SQUID rings. 

The correspondence principle in quantum mechanics is usually expressed in the form: "For those 
quantum systems with a classical analogue, as Planck 's constant becomes vanishingly small the expectation 
values of observables behave like their classical counterparts" [12]. for SQUID rings such an expression 
turns out to be problematic and we find that an alternative expression is more appropriate [sic] : "Consider 
h fixed (it is) and scale the Hamiltonian so that when compared with the minimum area h/2 in phase 
space: (a) the relative motion of the expectation values of the observable become large and (b) the state 
vector is localised. Then, under these circumstances, expectation values of operators will behave like their 
classical counterparts" [2]. 

In order to achieve localisation and model a dissipative chaotic-like system in its correspondence limit 
we need to introduce decoherence in the right way. Quantum state diffusion has proved particularly 
successful in many studies of non-linear system. Here we have an Ito increment equation for the state 
vector of the form IIOIIIII 
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dt + Y, 
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here Lj = ^/T^dj, where aj is the annihilation operator and dt and the d^ are complex Weiner increments 
satisfying d^"^ = = and d^d^,* = dt [lOtlllj where the over-bar denotes the average over infinitely 
many stochastic processes. 

QSD, however, is not the only unravelling of the master equation and for our results to be general they 
should be demonstrated to be independent of this choice. This point may be emphasised by observing 
that previous studies have shown that entanglement can be dependent upon the choice unravelling |13j . 
We therefore now choose another unravelling against which we may check our results. We choose an 
unravelling that is very different from QSD as it is based on a discontinuous photon counting measurement 
process - rather than a continuous interaction - namely quantum jumps jl4|,ll5j. Again this model takes 
the form of a stochastic Ito increment equation for the state vector but now of the form 
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where dNj is a Poissonian noise process such that dNjdN/^ = SjkdNj, dNjdt = and dNj 
i.e. jumps occur randomly at a rate that is determined by (.L^jLj 



L]LAdt 
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Figure 1: Mean entropy of entanglement as a function of (3 for the chaotic- like and periodic (entrained) 
states. Here we see that the entropy of entanglement for the system in the chaotic state does not vanish 
as /3 approaches the classical regime. (Note: Figure and caption reproduced from [1]) 

In [1] we studied the entanglement dynamics (characterised via the entropy of entanglement S {pi) = 
— Tr[pj In Pi]) in two coupled Duffing oscillators [T] (extending one dimensional analysis in, for example, [5l 
[8]). The Hamiltonian for each oscillator was given by 

Hi = ]^pi + ~ \li + ^ + Y ^'i'P' +Pi1i) (3) 

where qi and pi, Li = \/2Tiai (for i = 1,2), where is the annihilation operator. Here gi = 0.3 
and Fj = 0.125, [HEllH]. In this work the parameter /3 is a scaling parameter used to generate the 
correspondence limit. The Hamiltonian for the coupled system is: 

H = Hi+H2 + ^gi<?2 (4) 

with n = 0.2. 

The dynamics of the oscillators have two distinct modes of operation; entrained & periodic and 
un-entrained &: chaotic. When the oscillators are entrained we found that, as one would expect, the 
entanglement falls as the system approaches the classical regime. In the un-entrained &: chaotic mode of 
operation we found that significant average entangled was manifest both in the quantum and classical 
limit. These results are shown in Fig. [1] using quantum state diffusion and Fig. [2] for quantum jumps 
unravellings of the master equation. 

In [2l[3] we extended this investigation to SQUID's. Here the "classical" dynamics are described by 
the resistively shunted junction (RSJ) model: 

+ -- + ^ + I. sm j = /, sm (..t) (5) 

where $ is the magnetic fiux contained within the ring C, Ic: L, R, I^, ojd and $o = h/2e are the 
external flux bias, capacitance and critical current of the weak link, ring inductance, resistance, drive 
amplitude, drive frequency and fiux quantum, respectively. Here, C = 1 x 10~^^F, L = 3 X 10~^'^H, 
R = lOOf], /3 = 2, = ojQ, = 0.5^>o and h = 0.9 ^tA. 
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Figure 2: The calculation of figure [T] using quantum jumps instead of quantum state diffusion. Again we 
show the mean entropy of entanglement as a function of /3 for the chaotic-like and periodic (entrained) 
states. As with quantum state diffusion we see that when using quantum jumps the entropy of entan- 
glement for the system in the chaotic state does not vanish as /3 approaches the classical regime. (Note: 
Figure and caption reproduced from [l]) 



[sic [3]] "We can then rewrite ([5]) in the standard, universal oscillator like, form by making the 
following definitions: ujq = 1/^/LC, t = ujQt, ip = — ^^)/^q, ipx = ^x/^O: /3 = 27rL/c/*J>05 ^ = ^d/^o-, 
'Pd = IdL/^0 and ^ = l/2u;()RC. This yields the following equation of motion: 



(fip dip /3 

^ + 2C- + c,+ -sin[2vr(v, + ,,, 



Pd sm [UJT) 



(6) 



In this system of units we then see that we can scale the system Hamiltonian through changing either 
C —7- aC or L — 7- 6L so long as we also make the following changes: R — >• ^JhjaTi^ Id — Id/Vb and 
ujd uJdl^foh. . . .We change a so that C varies between 1 x 10~^^F (quantum limit) and 1 x 10~^ F 
(classical limit), changing other circuit parameters in line with the above methodology." 
The Hamiltonian is: 
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(7) 



with 



ih. 



As usual we define: Xi = sj Cujo/h^i and pi = ^Jl/hCoJoQi. and II[ = Hi/hwo so that 



fr/ _ Pi , [Xi 



cos (flXj) 



(8) 



where 



1/2 



One further correction to the Hamiltonian is needed to correctly introduce damping [5] which now 
becomes: 



2 



(9) 
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Figure 3: Mean entanglement entropy as a function of Capacitance two coupled SQUID rings using (a) 
quantum state diffusion and (b) quantum jumps unravellings of the master equation. In both figures we 
see that the entanglement entropy for system does not vanish even as it approaches its classical limit. 
Note: that unlike in Fig. [T] and Fig. [2] in this figure the quantum limit is on the left hand side and the 
classical limit on the right. 

So, for two coupled SQUID's we have 



where we have chosen fj, = 0.2 (as this is the value that we used in [T]). 

In Fig. O^a) we show the mean entanglement of the two SQUID rings found by using the Quantum 
state diffusion unravelling of the master equation (these results were also presented in [3]). Here small 
capacitance is the quantum limit and large capacitance is the correspondence limit. The capacitance was 
changed via use of the scaling parameters a of the discussion above, [sic [3]] "However we note that the 
entanglement entropies presented here are is the average entanglement over either a long time period or 
many similar trajectories. It is not the entanglement associated with the average density operator taken 
of many experiments. This average entanglement cannot therefore be considered usable in a quantum 
information sense. In figure we show this average entanglement entropy. Here the averaging of each 
trajectory was determined on a point by point basis. A sufficient averaging was used so as to ensure that 
the results presented here had settled to within a percent or so ... As for the Duffing oscillators, here the 
mean entanglement does not appear to vanish in the classical limit (large capacitance). Another surprising 
feature in common with the Duffing oscillator results is that the average entropy is not maximum at the 
most quantum limit (smallest capacitance)." 

In Fig. we present the result of this paper - here we have simply reproduced the calculations of 
Fig. El^a) using the quantum jumps unravelling of the master equation. We note that for the quantum 
jumps model that - especially in the quantum limit - it takes much longer for the averages to settle to 
their final values and there is some small error attached to each of the data points. However there is a 
good qualitative agreement between these results and those obtained for the Duffing oscillator. Is seems 
then that such persistent entanglement is independent of measurement process and that our results may 
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well be quite general in nature. 
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